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ABSTRACT 



A ten-level primitive equation ocean circulation model 
is used to investigate the formation and evolution of 
large-scale thermal anomalies observed in the central North 
Pacific Ocean during the fall and winter of 1976. Several 
initial value model integrations of 4 months duration are 
carried out in order to help explain the observed anomaly 
development. Initial and verifying ocean data down to 400 m 
depth are obtained from the NORPAX ships of opportunity pro- 
gram. Anomalous atmospheric wind forcing is obtained from 
Namias * monthly mean sea-level pressure anomalies, while 
climatological heat fluxes are used. 

The skill with which the model simulates the observed 
anomaly evolution in the different experiments is estimated 
synoptically and measured statistically by calculating root 
mean square (RMS) temperature errors and SI skill scores. 
Analysis indicates that anomalous atmospheric wind forcing 
improves the model predictions in the upper levels. For 
this particular winter case using climatological heating, 
however, knowledge of the initial anomalous temperature con- 
ditions does not improve the model results. The model skill 
at upper levels exceeds both persistence and climatology 
(forecast of zero anomaly) while at the lower levels it is 
comparable to persistence but not as good as climatology. 
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I 



INTRODUCTION 



Early investigations by Namias and others indicated that 
the winds are an important source of kinetic energy for the 
oceanic circulation. But just as the mean atmospheric cir- 
culation helps drive the mean large scale ocean circulation, 
anomalous atmospheric forcing can be assumed to drive anoma- 
lous currents in the ocean (Malkus, 1962). This premise was, 
in part, the motivating force for the present model study. 

The effects of anomalous atmospheric forcing on the oceans 
have been a matter of intuitive guessing for years, but few 
studies have attempted to assess in more than a descriptive 
or qualitative manner the nature or magnitude of these 
effects. A notable exception is the recent work of Davis 
(1976, 1977) . 

One of the most important manifestations of the effects 
of anomalous atmospheric forcing on the oceans is the devia- 
tion of the ocean thermal structure from the long term cli- 
matological norm. The space and time scales of these devia- 
tions or anomalies are known in a very general sense, but 
the exact generating mechanism or mechanisms remain unclear 
at present. It is the intent of this investigation, within 
the larger context of the Anomaly Dynamics Study (ADS), to 
"Describe and understand the mechanisms responsible for 
large scale climatological anomalies and related changes in 
the thermal structure and the general circulation of the 
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upper layers of the ocean in order to predict their forma- 
tion and evolution.” (Kirwan, 1976) The contribution of 
anomalous atmospheric wind forcing as a generating mechanism 
for changes in the thermal structure of the ocean will be 
emphasized . 

Previous model investigations of the dynamics of large 
scale temperature anomaly formation in the ocean have been 
limited by the lack of verifying data, at the surface and 
at depth, over sufficiently large areas. In cases where 
anomalous atmospheric forcing for a certain period of time 
was added to climatological wind forcing (Namias, 1965; 
Arthur, 1966; Jacob, 1967), model results show an important 
influence on the ocean by the anomalous forcing. 

In the present study, an ocean circulation model is 
used to make a number of initial value predictions of the 
evolution of large scale ocean anomalies.. The different 
model predictions are designed to examine the effect on the 
prediction of including or excluding certain physical pro- 
cesses in the model. The study is unique in that a care- 
fully collected and analyzed set of ocean temperature data 
are used. 

It must be emphasized at this point that this study is 
not meant to be a test of the predictive capability of the 
model itself. First, a very important part of the atmos- 
pheric forcing fields, namely anomalous heat flux that is 
required for the testing of the prediction capabilities of 
any ocean model is not available on the time or space scales 
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required. Consequently, only anomalous winds are used in the 
present studies. Second, ocean modeling itself is still 
relatively unsophisticated and certainly untested in model- 
ing large scale thermal anomalies observed in the oceans 
even if all the atmospheric forcing were known. In addition, 
this study will serve as an initial investigation into the 
adequacy of currently available temperature data in describ- 
ing the true thermal condition of the Central Pacific. With 
the current oceanic and atmospheric data available and with 
the current state of the science of ocean modeling, all that 
can be tested are a number of hypotheses regarding the rela- 
tive importance of various dynamic forces that act to drive 
the oceans. Any conclusions resulting from this study then 
are subject to re-evaluation as ocean modeling and monitoring 
capabilities improve. The 1976-77 fall-winter period was 
selected for this study because of the large amplitude anoma- 
lies in the atmosphere and ocean that existed in the North 
Pacific at that time. 
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II 



DATA PROCESSING 



The modeling experiments described below made use of 
anomalous atmospheric winds and oceanic temperatures. The 
first step in preparing the atmospheric data for input into 
the ocean model was calculation of sea-level pressure anoma- 
lies or departures from normal for all of 1976 and January 
and February of 1977 . All the atmospheric data were generous- 
ly made available to us by Dr. J. Namias of Scripps. Long 
term (20 year) mean values of sea level pressure for each of 
these months, as determined by Namias, were subtracted from 
monthly means in order to obtain the sea level pressure de- 
partures from norm. Analysis of these anomaly fields indi- 
cated the period from September 1976 until January 1977 ex- 
hibited sufficiently large anomalous sea level pressures to 
be suitable for our investigation. 

In order then to input Namias* anomalous pressure data 
into the model, it was necessary to extrapolate and then 
interpolate the values from Namias* 5®* latitude by 5® longi- 
tude (9x19) fields encompassing 20®N to 60®N and 145®E to 
125®W to the 33x33 fields extending from the equator to 65®N 
and 145®E to 125‘^W required by the model (see figure 1). 

The extrapolation of the anomalous pressure data toward the 
northern boundary assumed a slope that decreased by 50 per- 
cent at each grid point. Extrapolation towards the equator 
was such that anomalous u and v components of the wind field. 
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calculated from the anomalous pressures, go to zero at the 
equator. Once the extrapolation was completed, data north 
of 50^ N was smoothed several times in the east-west direc- 
tion in order to remove any small scale features introduced 
by the extrapolation. This procedure was carried out for 
the months of September 1976 through January 1977. 

With the 33x33 field of sea level pressure anomalies 
computed, the corresponding anomalous geostrophic u and v 
components of the wind were calculated. These anomalous geo- 
strophic wind fields were the only anomalous atmospheric 
conditions used in the ocean model. These fields (figure 2) 
are described in the Results section below. 

The next step in data preparation for input into the 
model was to compute the temperature anomalies that would 
serve as initial and verifying conditions for the model. 

All temperature data was provided by TRANSPAC, the ships of 
opportunity cooperative research effort. The TRANSPAC ex- 
pendable bathythermograph (XBT) field program is a part of 
the North Pacific Experiment (NORPAX) which is charged with 
investigating how oceanic variability influences the weather 
and climate, and in turn how the weather and climate varia- 
bility influence the oceanic circula tions . Prior to TRANSPAC, 
data on the thermal structure of the ocean was inadequate, 
thus leading to the inception of the TRANSPAC XBT field pro- 
gram to provide the subsurface temperature data necessary 
for NORPAX. The observational phase of the program began in 
January 1975 utilizing ships of opportunity (i.e. merchant 
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ships voluntarily providing a no-cost observation platform) 
that routinely traverse the sea lanes between the United 
States and Japan. The teraperature anomaly data provided 
through TRANSPAC was obtained by subtracting climatology 
means from the monthly mean observed temperatures. Anomalous 
temperatures were provided at ten levels including 0,30, 60, 

90, 120, 150, 200, 250, 300 and 400 meters. For initializ- 

ing temperature fields in the model, the anomalous tempera- 
ture field for September 1976 was added to the model gener- 
ated climatology for September. Calculation of verification 
statistics entailed use of anomalous temperature fields for 
October through December of 1976 and January 1977. Again, 
the TRANSPAC anomaly data had to be processed for use in the 
model because the domain and standard levels of the TRANSPAC 
data- did not coincide with the domain and standard levels of 
the model. It was decided, in order to avoid major changes 
to the model, to extrapolate from the 2® latitude by 5° longi- 
tude (11x15) TRANSPAC fields extending from 30°N to 50°N 
and 160°E to 130°W to the area and gridpoints encompassed 
by the model (see figure 1) . 

The first step in preparing the TRANSPAC temperature anom- 
aly data was to interpolate the data from the ten levels of 
TRANSPAC to 11 intermediate levels of the model, specifically 
to 0, 20, 45, 80, 125, 200, 325, 600, 1200, 2200 and 4000 

meters. Below 400 meters, values were assigned to the inter- 
mediate levels by decreasing the 400 meter anomalous tempera- 
ture to a value of zero at the 2200 meter intermediate level. 
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This was accomplished by multiplying the 400- meter value by 
a cosine function which decreased with depth* For example, 
if T*(400) is the given TRANSPAC anomaly at any grid point, 
then the anomaly at depth Z, T' (Z) , was specified by 

T'(Z) = T'(400) cos [( 3 ^^) . 

This is a very important assumption, since it assumes the 
temperature anomaly cannot change sign with depth below 
400 meters. The resulting temperatures at the 11 intermedi- 
ate levels were then averaged between any two levels to de- 
fine the temperature at the ten prognostic levels in the 
model which are located at 10.0, 32.5, 62.5, 102.5, 162.5, 

262.5, 462.5, 900.0, 1700.0 and 3100.0 meters. 

The next step in preparing the TRANSPAC data was to ex- 
pand it to cover the entire area of the model. This was 
done by expanding the temperature anomaly at the border of 
the 11x15 TRANSPAC area outward to zero at the outer boundary 
of the model. In order to provide a smooth transition into 
the area of no data, the extrapolation was along a sine func- 
tion as shown in figure 1. Areas of the array which could 
be filled by both longitudinal and latitudinal extrapolation 
were assigned a temperature anomaly value which was an 
average of the two extrapolations. Once the entire geographi- 
cal area of the model was filled on a 2 ® latitude by 5® 
longitude grid, an interpolation scheme was used to fill 
the 33x33 array of model grid points. The 33x33 array of 
anomalous temperatures at the ten levels of the model for 
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September 1976 was then added to the model generated clima- 



tology for Septemb 
ture conditions of 
October 1976 to Ja 
for computing veri 
mode 1 integrations 



er 15 to constitute 
the model. Similar 
nuary 1977 served as 
fication statistics 
as described below. 



the initial tempera- 
33x33x10 arrays for 
the verifying data 
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III. MODELING APPROACH 



In order to test a variety of atmospheric and oceanic 
processes that may be responsible for the formation and 
evolution of the large-scale temperature anomalies that are 
of interest, a series of initial value model integrations 
have been carried out using a 10-level primitive equation 
model in a closed rectangular basin. The model is that of 
Haney (1974) which has been improved to include time de- 
pendent seasonal forcing by the atmosphere, a parameteriza- 
tion of surface wind and convective mixing (Haney and 
Davies, 1976), and nonlinear lateral eddy viscosity based 
on two-dimensional turbulence theory (Haney and Wright, 
1975). The ocean model is driven by wind and differential 
heating which is calculated from the model sea-surface 
temperature and prescribed values of the atmospheric solar 
radiation, total cloudiness, surface air temperature, rela- 
tive humidity and winds. These quantities are made up of 
a climatological part and in some cases an anomalous part. 

To initialize the model for a prediction experiment, 
the model is first integrated over a long ”spin-up’* period 
(240 years) using climatological atmospheric forcing. The 
last year of this integration is referred to as the "model 
climatology" and the resulting currents and temperatures, 
especially in the mid-ocean areas, bear a close resemblance 
to those in the North Pacific Ocean. The initial conditions 
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for the prediction experiments are then obtained by combining 



a given 


anoma 


model c 


limato 


time of 


the y 


One 


e the 


anomaly 


field 


runs, e 


ach wi 


As stat 


ed in 


model e 


xperim 


various 


a tmos 


and evo 


lu tion 


The 


f ive 


scribed 


by wr 



thermodynamic 
Table I) , 




dt p 

o 




ly, whether oceani 
logy for the same 
ear . 

anomalous atmosphe 
s were prepared fo 
th different model 
the introduction, 
ents was to examin 
pheric and oceanic 
of the observed a 
different model in 
iting the equation 
energy equation a 

V„p - fK X + V* 
, T) + V* (aVt) + -^( 



c or atmospheric, with the 
parameter for the appropriate 



ric forcing and temperature 
r input into the model, five 
configurations, were made, 
the purpose of the different 
e the relative importance of 
processes in the formation 
nomalies . 

tegrations are most easily de~ 
of horizontal motion and the 
s follows, (symbols defined in 



(AW^) 



+ 




9V 

Tz 



H 



( 1 ) 



K 



il) 



( 2 ) 



The terms on the right hand side of (1) are the pressure 
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of surface wind and convection, horizontal diffusion and 
vertical diffusion. The ocean model is driven by a downward 
flux of heat and momentum (stress) at the surface which 
enters through boundary conditions in the vertical diffusion 
terms : 
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8 Z^Z = 0 = Q = downward heat flux at the surface 

8V 

H 

(K - 5 -=— ) _ = T = surface stress 

d ^ Z = U 

In all of the model experiments climatological values of the 
heat flux, 0 , were used. The stress/ i / was calculated from 

''s ' 

where is a prescribed surface wind which is the sum of 

a climatological part and an anomalous part; the latter 

being calculated from Namias * sea level pressure anomalies 

described above. Thus, only V has an anomalous part, say 

V * and it enters the model equations in two places; as a 
o 

boundary condition (v/ind stress) in the vertical diffusion 

term, (K . ; and directly in the vertical mixing term 

d Z Z— 0 

M(Q,T). The ability exists in the model then to allow the 
anomalous surface winds to enter into either, both or 
neither of the terms. In this way, we examine separately 
the role of anomalous wind driven currents and . anomalous 
wind generated vertical mixing in the generation of oceanic 
thermal anomalies . 

In physical terms, the first model configuration, model 

1 , describes the free evolution of an initial temperature 

anomaly since V ' = 0 in this model. In model 2, the 

S 

effect of the anomalous wind stress in driving anomalous 
Ekman type flow is added, since Vg ' = 0 only in the mixing 
term in this model. According to steady Ekman theory this 
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flow would be 90° to the right of the anomalous surface 
winds. Model 3 incorporates the effects of anomalous wind 
stress in both Ekman type flow and surface wind mixing. 

This mixing process is expected to produce increased verti- 
cal mixing in regions of strong anomalous winds. Model 4 
is similar in its dynamics to model 2 but has the initial 
anomalous temperature fields equal to zero. Likewise, 
model 5 is similar to model 3 but with the initial anomalous 
temperatures equal to zero. These last two experiments are 
designed to evaluate the importance of having good initial 
conditions when making a prediction. 

In order to provide a quantitative comparison between 
the different models, a variety of statistics were computed 
for the five model cases at 30-day intervals of time. The 
statistics calculated included the average absolute predic- 
tion error, root mean square error and an SI score associated 
with the model generated temperature anomalies. The follow- 
ing definitions indicate exactly what was measured by the 
statistics. Average absolute error and root mean square 
error are defined as follows: if in a series of N forecasts 

(i.e. N different gridpoints), represents the i-th fore- 

cast and 0^ the corresponding observation, the average 
absolute error is given by 



e = 



Z F. - 0. 

* 1 ^ 

N 
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and the root mean square error by 



RMSE 





2 



(Brier and Allen, 1951) . 

The SI score is a dimensionless measure of the difference 
between predicted and observed horizontal gradients and is 
given by 



(Teweles and Wobus, 1954), 

where e^ is the error in the forecast temperature anomaly 
gradient in the horizontal and G is the observed tempera- 

li 

ture anomaly gradient. The SI score was further divided in 
order to measure the skill in forecasting both the north- 
south and the east-west gradient of the temperature anomaly. 

The various model predicted and TRANSPAC derived tem- 
perature anomaly fields were compared beginning with the 
October 1976 fields at 30-day intervals of time and at nine 
levels of the model. Statistics were computed only over an 
area from 35®N to 45®N and 130°W to 180®. The verification 
was restricted to the region east of the dateline because 
the region west of the dateline is generally characterized 
by numerous mesoscale temperature features (Bernstein and 
White, 1977). Since the intent of this investigation was in 
part to study present capabilities in modeling large scale 
thermal anomalies, it was assumed the model could not effec- 
tively predict such mesoscale features. 



21 



IV 



RESULTS AIID CONCLUSIONS 



A vast amount of data is made available by this study 
and as is always the case, interpreting such a volume of 
data in light of the original intentions of the study is 
difficult. Describing the actual evolution of the thermal 
anomalies in the observed and various model cases is perhaps 
the most logical starting point. A later more detailed 
treatment of the verification statistics themselves will 
serve to highlight various aspects of the described results. 

Since the effects of anomalous wind stress in producing 
large scale thermal anomalies is considered paramount, a 
description of the anomalous atmospheric forcing for 
September 1976 to January 1977 should precede any descrip- 
tion of observed or model evolutions of thermal anomalies. 
Figure 2 shows the anomalous sea-level pressure fields for 
these months, contoured at 2 mb intervals with shaded areas 
indicating negative anomalous pressures. These anomalous 
pressure fields are shown only over the area of the tem- 
perature fields given in figures 3 through 8. The September 
1976 mean circulation over the central Pacific was dominated 
by a trough at 700 mb which had progressed to a position 
south of the Aleutians. Downstream, the western North 
American ridge became established as an Arctic ridge 
(Taubensee, 1976). This mean 700 mb trough in the central 
Pacific is evidenced at the surface by negative anomalous 
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pressure. In October 1976/ the mean 700 mb trough over the 
eastern Pacific retrograded somewhat, as did the downstream 
ridge over western North America (Wagner, 1977a). The nega- 
tive anomalous pressure pattern at sea level is seen in 
figure 2 to move westward in October, being replaced to the 
east by positive anomalous pressure associated with the 
ridge over western North America. The mean 700 mb trough 
over the central Pacific remained relatively stationary but 
intensified significantly (2 standard deviations from the 
mean) from November 1976 to January 1977 (Dickson, 1977; 
Taubensee, 1977; Wagner, 1977b). This intensification is mani- 
fested in the sea level pressure by an extreme deepening 
of the negative anomalous pressure fields for November 
through January. 

The observed anomalous temperature data for September 
1976, November 1976 and January 1977 at 10, 102, 162 and 262 

meters are shown in figure 3. The observed evolution during 
this winter period in the upper level is characterized by 
extreme cooling over all but the easternmost edge of the 
area. In the mid-levels (102 and 162 meters), the September 
field shows an extremely intense disturbance in the monthly 
mean temperature anomaly at and just east of the dateline. 

This disturbance propagates westward out of the area by 
November. The 102 meter level exhibits a general cooling 
over the winter months, similar to that in the upper level. 

At 162 meters there is a general trend towards warming with 
the axis of temperature anomalies shifting from generally 
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north-south to east-west. The 262 meter field shows a less 
intense anomaly structure than the mid-levels for September 
but exhibits a more intense warming along the east-west axis 
than had occurred by January at 162 meters. 

Figure 4 shows the anomaly evolution predicted by model 1 
(climatological winds). The dominant feature in the upper 
levels is the development of an intense, basically mesoscale 
feature in the western portion of the area. This development 
is due to the strong north-south advection of mean tempera- 
ture by the anomalous geostrophic surface currents induced 
in the model by the intense disturbance in the western part 
of the basin at mid-levels. This same type of development 
occurs in the middle and lower levels so that a similar but 
less intense structure develops in the vertical. These 
features tend toward a slow westward propagation during the 
winter months, indicating partial development and movement 



as free Rossby waves. 








Model 2 


introduces 


anomalous 


winds, but only in 


the sur- 


face stress 


term (clima 


tological 


winds are used in 


the verti 



cal mixing term). Including this effect produces anomalous 
Ekman type flow in the upper levels which definitely improves 
the model predictions by cooling the upper layers (see figure 
5). The intense mesoscale features in the west produced by 
the strong north-south advection of mean temperature de- 
scribed in model 1 are still very apparent in the anomaly 
structure however. These anomalously cold mesoscale features 
are intensified and expanded in their horizontal scale in the 
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upper levels by the superposition of the large scale cooling 
resulting from the anomalous wind stress. In the middle to 
lower levels there is a general warming, similar to the ob- 
served case, which tends to reduce the intensity of the cold 
eddies . 

Model 3 includes the effects of anomalous winds in both 
the surface stress term and the vertical wind mixing of heat. 
This wind mixing process would be expected to produce in- 
creased vertical mixing in regions of increased surface winds 
but as seen in figure 6, the results for model 3 do not vary 
significantly from model 2. In the ocean model, surface wind 
mixing depends in part on the net surface heat flux. If the 
net flux is downward, wind mixing occurs, but if the net sur- 
face heat flux is upward as it is in this winter case there 
is little or no wind mixing. Model 3 results then show only 
slight differences from model 2 in the uppermost layers. 

Model 4 includes anomalous winds only in the surface 
stress term (like model 2) but the initial temperature anomaly 
is set to zero in order to evaluate the importance of knowing 
the initial conditions. In this particular winter case of 
very strong anomalous wind forcing, the temperature anomalies 
in the upper levels are predicted slightly better than in 
models 2 or 3, with general cooling over most of the area as 
seen in figure 7. Results in middle to lower levels are 
worse, however, with the 102 m level showing development of 
small scale warm anomalies near 30-35N where cooling actually 
occurred and the 262 m level showing cooling to the northeast 
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where warming is observed. It may be noted at this point 
that the superposition of the results for model 4 with 
initial anomalous temperature equal zero and model 1 with 
no anomalous wind forcing will give almost exactly the 
results observed in model 2. This is cldar evidence that 
in the model, the effects of external (wind) forcing and the 
effects of internal dynamical development take place rather 
independent of one another. 

Model 5 is the same in its dynamics as model 3 but again 
the initial anomalous temperature is set to zero. As pre- 
viously discussed, models 2 and 3 show little variation and 
as seen in figure 8, model 5 results vary only slightly from 
model 4. 

The actual verification statistics, of course, provide more 
concrete evidence to the similarity and differences between 
the various models in successfully describing the evolution 
of the large scale temperature anomalies as observed. Table 
n presents these statistics in the following format: nine 

levels from 10 meters to 1700 meters depth are considered; 
the ratio of the root mean square error (RMSE) to the RMSE 
of ' climatology (a prediction of zero temperature anomaly 
everywhere) is given as the first number followed by a 
similarly non-dimens ional ized SI score for 30 day intervals 
from October 1976 to January 1977. A ratio of RMSE * s and Sl*s 
were used as a measure of skill for easy interpretation since 
a ratio less than one indicates skill over climatology 
whereas a ratio greater than one indicates less skill than 



26 



climatology. The skill of a forecast of persistence is also 
shown. These ratios for 10, 102, 162, and 262 meters have 

been plotted as a function of time in figures 9a and 9b for 
easy reference. 

The plot of skill as a function of time at 10 meters 
seen in figure 9a indicates an extremely rapid loss of skill 
as early as 30 days, followed by improved skill in the models 
(2, 3, 4 and 5) having wind forcing. Model 1 with no 

anomalous wind forcing shows no skill after 30 days but 
models 2 and 3 clearly demonstrate the benefits of consider- 
ing such anomalous forcing. Surprisingly, models 4 and 5 
with zero initial anomalous temperature conditions show con- 
stant improvement with time even over models 2 and 3. This 
would apparently indicate the greater importance of knowing 
the anomalous wind forcing than of knowing the initial tem- 
perature field exactly, especially in a case as this one 
where the anomalous winds were very strong. 

Perhaps the most conspicuous result is the fact that the 
model predictions at the surface are better when the initial 
anomaly is zero (model 4 and 5) than when it is the observed 
anomaly (model 2 and 3). In the latter case, strong geo- 
strophic flow in the north-south direction is induced at 
the surface and this flow advects the mean isotherms in the 
model (cold to the north and warm to the south) to produce 
a large amplitude disturbance in the western third of the 
basin which does not exist in the observed pattern. As de- 
scribed earlier, the strong meridional flow is due to the 
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intense anomaly in the initial temperature field at and be- 



low 102 meters* One possible explanation for the model 
discrepancy is that the model simply does not propagate this 
developing disturbance westward and out of the domain fast 
enough. However, an examination of the surface anomaly 
data west of the dateline shows no large disturbance there 
to which the predicted disturbance in figure 4 could corre- 
spond if the model error was simply one of phase error. 

A second possibility is that in September there existed a 
compensating deep temperature anomaly (below 400 meters) 
whose induced vertical shear nearly cancels that induced by 
the anomaly observed between 102 and 262 meters. In this 
case, the model assumption that the initial anomaly decreases 
slowly to zero at depth, without changing sign (see DATA 
PROCESSING) is wrong. If this second explanation is true, 
it implies that in order to predict or understand sea surface 
temperature anomalies, one must measure the anomaly field 
all the way down to the bottom of the sea. Unfortunately, 
it will require many case studies of this type before this 
question is adequately resolved. 

The plot of skill versus time for 102 meters in figure 
9a shows a somewhat less rapid loss of skill initially but 
demonstrates a fairly high amplitude oscillation about 1.0 
in the 30 to 120 day period of time. Comparison of the 
curves for models 4 and 5 and models 1, 2 and 3 indicates 

a relation between these oscillations and the initial tem- 
perature conditions since this alternating behavior is not 
displayed by models 4 and 5. 
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Figure 9b shows an extreme and rapid loss of skill (com- 
pared to climatology i.e. a forecast of zero anomaly) for 
persistence and the models at 162 and 262 meters. An oscil- 
lation in performance characterizes the 162 meter graph. 

It is important to note also that while models 4 and 5 have 
greater skill, they are not as good as climatology. The 262 
meter plot does not exhibit this oscillatory characteristic, 
presumably as a result of less intense anomalous initial 
temperatures at this depth. It is rather disheartening that 
the model loses skill so fast at these levels. It is sur- 
prising too that persistence is also a very bad forecast 
at these depths which are generally believed to undergo 
slow changes. 

The SI scores computed to measure skill in predicting the 
gradient of large scale thermal anomalies are shown in Table 
I. Models 4 and 5 consistently demonstrate greater skill 
than either persistence or the other models in predicting 
gradients. This indicates of course that the ocean model has 
little skill in computing gradients since any field more 
closely resembling a zero field possesses greater skill. If 
the model SI scores are broken up into north-south and east- 
west components, however, there is greater relative skill in 
predicting north-south gradients than east-west gradients. 

The results of this study can be summarized as follows. 

In the upper levels, the models which include anomalous wind 
forcing show skill over climatology out to 4 months, i.e. the 
ratio of RMSE of the model to RMSE of climatology is less 
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than one. However, knowledge of the initial anomalous tem- 
perature conditions degrades the prediction skill. Two 
possible explanations were offered but a variety of addi- 
tional case studies are needed to properly resolve this 
problem. In the lower layers (162 to 262 meters), none of 
the models show skill over climatology and consideration 
of neither anomalous wind forcing nor initial anomalous tem- 
peratures improves the models predictive capabilities. The 
rather poor model results at these depths are a serious con- 
cern because it is not likely that the addition of anomalous 
heat fluxes to the model predictions will have much effect 
at these levels. In the future, it is planned to repeat 
this type of study for cases in which the anomalous atmo- 
spheric forcing was very weak and also to repeat the present 
study adding anomalous heat fluxes being computed by other 
NORPAX investigations . 
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10 M (RMSE RATIO/Sl RATIO) 


OCT NOV DEC JAN 




PERS 




1.17/1.53 


1.11/1.42 




1.12/1.61 




1.18/1.76 




Ml 




1.12/1.46 


1.08/1.61 




1.07/1.77 




1.11/1.79 




M2 




1.04/1.48 


0.89/1.58 




0.70/1.56 




0.77/1.52 


-- 


M3 




1.06/1.53 


0.87/1.66 




0.68/1.57 




0 .77/1.51 


M4 




0 .92/0.99 


0 .82/1.05 




0.63/0.99 




0.69/1.01 




M5 




0.96/1.13 


0.80/1.14 




0.62/0.99 




0 . 69/0 . 99 





32 M (RMSE RATIO/Sl RATIO) 


OCT NOV DEC JAN 


PERS 




1.18/1.53 




1.12/1.50 




1.12/1.47 




1.13/1.64 




HI 




1.14/1.48 




1.06/1.68 




1.07/1.77 




1.08/1.85 




M2 




1.16/1.63 




1.01/1.82 




0.81/1.56 




0.83/1.55 




M3 




1.08/1.91 




0.98/1.90 




0.79/1.57 




0 . 83/1 . 54 




M4 




1.02/1.08 




0 .96/1. 38 




0.75/0.99 




0.77/0.95 




M5 




0.93/1.43 




0.93/1.45 




0.73/0.98 




0.77/0.93 





62 M (RMSE RATIO/Sl RATIO) 


OCT NOV DEC JAN 


PERS 




1.35/1.23 




0.97/1.80 




0 .94/1.43 




0.99/1.83 




Ml 




0.96/1.15 




0.99/1.70 




0.98/1.56 




1.08/1.91 




M2 




0.99/1.21 




1.05/1.94 




0.94/1.75 




0.81/1 .79 




M3 




1.04/1.36 




1 . 07/2 . 01 




0.95/1.74 




0.81/1.79 




M4 




1.01/0.97 




1.05/1.29 




0.98/1.36 




0.74/1.20 




M5 




1.06/1.09 




1.07/1.41 




0.99/1.36 




0.74/1.21 





Table II 
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102 M (RMSE RATIO/Sl RATIO) 


OCT NOV DEC JAN 




PERS 




0.87/1.22 




1. 38/1.67 




0 . 84/1.63 




0.95/2.29 




Ml 




0.91/1.17 




1.14/1.52 




0.79/1.32 




1.03/1.94 




M2 




0.91/1.18 




1.17/1.55 




0.86/1. 45 




0 . 98/2.01 




M3 




0.94/1.24 




1.24/1.58 




0 . 91/1 . 50 




0.98/2 .02 


M4 




1.00/0.99 




1.03/1.05 




1 .05/1 . 21 




0 . 97/1.31 




M5 




1 .06/1 .10 




1 . 17/1 . 14 




1.09/1. 26 




0.98/1.32 





162 M (RMSE RATIO/Sl RATIO) 


OCT NOV DEC JAH 


PERS 




1.04/1.46 




1.85/2.07 




1 . 26/1 . 80 




1.88/1.77 




Ml 




1.00/1.31 




1.54/1.91 




0.91/1.34 




1.59/1.66 




M2 




1 . 00/1 .32 




1.46/1.89 




0.91/1.53 




1.65/1.74 




M3 




1 . 00/1. 32 




1.46/1.89 




0 .91/1 .53 




1.65/1.74 




M4 




1 . 00/1.03 




1 . 00/1 . 10 




1.00/1.15 




1.06/1 . 18 




M5 




1.00/1.03 




1.00/1.10 




1.00/1.15 




1.06/1.18 





262 M (RMSE RATIO/Sl RATIO) 


OCT NOV DEC JAN 


PERS 




1.64/1.76 




1.73/2.29 




1.69/2.10 




0.97/1.38 




Ml 




1 .57/2.03 




1.55/2.11 




1 .38/1.85 




1.03/1.38 




M2 




1. 50/2.02 




1.45/2.09 




1.31/2.00 




1.05/1 .45 




M3 




1. 50/2.02 




1 . 45/2 .09 




1. 31/2.00 




1.05/1.45 




M4 




1.00/1.02 




L. 00/1. 13 




1 . 00/1.18 




1.03/1.11 


M5 




1.00/1.02 




1.00/1.13 




1.00/1. 18 




1.03/1.11 





Table II (cont) 
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462 M (RMSE RATIO/Sl RATIO) 


OCT NOV DEC 


JAN 




PERS 




1.62/1.87 


1 .00/1 . 66 




1.06/1.81 




0.89/1.37 




Ml 




1.62/1.78 


0 . 88/1. 35 




0.94/1.38 




0.86/1.17 




M2 




1.54/1.77 


0 .88/1.36 




0.89/1.41 




0 .89/1 .24 




M3 




1.54/1.77 


0.88/1.36 




0.89/1.41 




0 .89/1 .24 




M4 




1.00/0 . 98 


1.00/1.03 




1.00/1.07 




1.06/1.11 




M5 




1.00/0 . 98 


1.00/1.03 




1.00/1.07 




1 . 06/1 .11 





900 M (RMSE RATIO/Sl RATIO) 


OCT NOV DEC JAN 


PERS 




1.45/1.68 




0.93/1.44 




1.00/1 .64 




0.86/1 .31 




Ml 




1.45/1.58 




0.79/1.19 




0.93/1.29 




0.86/1.11 




M2 




1.45/1.56 




0 . 79/1 . 18 




0.93/1.30 




0 . 86/1 . 16 




M3 




1.45/1.56 




0.79/1.18 




0.93/1. 30 




0.86/1.16 




M4 




1.00/0.99 




1.00/1. 02 




1 . 00/1 . 06 




1.04/1.08 




M5 




1 .00/0.99 




1.00/1.02 




1.07./1.06 




1.04/1.08 





1700 M (RMSE RATIO/Sl RATIO) 


OCT NOV DEC JAN 


PERS 




1.50/1.68 




0.80/1.44 




1.00/1.64 




0 .82/1 .31 




Ml 




1.50/1.58 




0.80/1.18 




0.83/1.29 




0 .82/1 .11 




M2 




1.50/1.57 




0.80/1.17 




0.83/1.30 




0.82/1.13 




M3 




1.50/1.57 




0.80/1.17 




0.83/1.30 




0 . 82/1.13 




M4 




1.00/0.99 




1.00/0.99 




1.00/0.99 




1.00/1.02 




M5 




1.00/0.99 




1 . 00/0 .99 




1 . 00/0 .99 




1 . 00/1 . 02 





Table II (cont) 



44 



I 



1 






LIST OF REFERENCES 



Arthur, R. A., "Estimation of Mean Monthly Anomalies of Sea- 
Surface Temperature," Journal of Geophysical Research , 

V. P- 2689-2690, 1966. 

Bernstein, R. L. and White, W. B., "Zonal Variability in the 
Distribution of Eddy Energy in the Mid-Latitude North 
Pacific Ocean," Journal of Physical Oceanography , v. 7_, 

P. 123-126, 1977. 

Brier, G. W. and Allen, R. A., "Verification of Weather Fore- 
casts," Compendium of Meteorology , p. 841-848, 1951. 

Davis, R. E., "Predictability of Sea Surface Temperature and 
Sea Level Pressure Anomalies over the North Pacific 
Ocean," Journal of Physical Oceanography , v. p. 249- 

266, May 1976. 

, "A Search for Short Range Climate Predictability," 

Journal of Physical Oceanography , (in press), 1977. 

Dickson, R. R., "Weather and Circulation of November 1976," 
Monthly Weather Review , v. 105 , n. 2, p. 239-244, 

February 1977. 

Haney, R. L., "A Numerical Study of the Response of an 

Idealized Ocean to Large-Scale Surface Heat and Momentum 
Flux," Journal of Physical Oceanography , v. n. 2, 

p. 145-167, April 1974. 

, and Davies, R. W., "The Role of Surface Mixing in 

the Seasonal Variation of the Ocean Thermal Structure," 
Journal of Physical Oceanography , v. n. 4, p. 504- 

510, July 1976. 

, and Wright, J., "The Relationship between the Grid 

Size and the Coefficient of Nonlinear Lateral Eddy 
Viscosity in Numerical Ocean Circulation Models," 

Journal of Computational Physics , v. J^, p. 257-266, 1975. 

Jacob, W. J., "Numerical Semiprediction of Monthly Mean Sea 
Surface Temperature," Journal of Geophysical Research , 

V. 72_, p. 1681-1689, 1967 . 

Kirwan, A. D., Memorandum dated July 9, 1976, Subject: 

Anomaly Dynamics Study Summary and Five Year Plan , 

Texas A and M. 

Malkus, J. S., "Interchange of Properties between Sea and 

Air," The Sea , v. Interscience Publishers, p. 88-294, 

1962. 



45 



E 




s 



Namias, J., "Macroscopic Association between Mean Monthly 
Sea-Surface Temperature and the Overlying Winds," 

Journal of Geophysical Research , v. 7_0, p. 2307-2318 , 

1965 . 

Taubensee, R* E., "Weather and Circulation of September 1976," 
Monthly Weather Review , v. 104 , n. 12, p. 1631-1637, 
December 1976. 

, "Weather and Circulation of December 1976," Monthly 

Weather Review , v. 105 , n. 3, p. 368-373, March 1977. 

Teweles, S. and Wobus , H. B., "Verification of Prognostic 

Charts, " Bulletin of the American Meteorological Society , 
V. 3_5, n. 10, p. 455-463 , December 1954 . 

Wagner, A. J., "Weather and Circulation of October 1977," 
Monthly Weather Review , v. 105 , n. 1, p. 553-560, 

January 1977a. 

, "Weather and Circulation of January 1976," Monthly 

Weather Review, v. 105 , n. 4, p. 121-127, April 1977b. 



46 



1 ^ 



I 



INITIAL DISTRIBUTION LIST 



1. Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 

2. Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93940 

3. Dr, G, J, Haltiner, Code 63Ha 
Department of Meteorology 
Naval Postgraduate School 
Monterey, California 93940 

4. Dr. R. L. Haney, Code 63Hy 
Department of Meteorology 
Naval Postgraduate School 
Monterey, California 93940 

5. LT. Wayne S. Shiver USN 
COMPHIBGRU ONE 
FPO San Francisco 96601 

6. Director, Naval Oceanography and Meteorology 1 
Building 200 

Washington Navy Yard 
Washington, D. C. 20374 

7. Meteorology Department 1 

Code 63 Library 

Naval Postgraduate School 
Monterey, California 93940 

8. LT. K. H. Hunt 1 

Fleet Weather Central 

McAdie Building 

NAS, Norfolk, Virginia 23502 

9. Dr. Joe Huang 1 

Great Lakes Environmental Research Lab-NOAA 
2300 Washtenaw Avenue 

Ann Arbor, Michigan 48104 

10. Prof. W. L. Gates, Chairman 1 

Department of Atmospheric Sciences 

Oregon State University 
Corvallis, Oregon 97331 

11. Dr. J. Namias 1 

Department of Oceanography 

Scripps Institution of Oceanography 
La Jolla, California 92093 



No. Copies 
2 

2 

1 

1 

2 



47 




1 



173ZQ5 



Thesi s 

CI 1775 Shiver . 

^ 1 Dynamical numerical 

prediction of large-^ 
scale thermal anomalies 
in the North Pacific 
Ocean, 



I 



Thesis 173205 

S4775 Shiver 

c.l Dynamical numerical 

prediction of large- 
scale thermal anomalies 
in the North Pacific 
Ocean. 




J 



thesS4775 

Dynamical numerical prediction of large- 




3 2768 001 95341 7 

DUDLEY KNOX LIBRARY 















3 







